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ABSTRACT 

We use complete enumeration and Monte Carlo techniques to study self-avoiding 
walks with random nearest-neighbor interactions described by voqiQj, where qi = ±1 is a 
quenched sequence of "charges" on the chain. For equal numbers of positive and negative 
charges (A^+ = A^-), the polymer with vq > undergoes a transition from self-avoiding 
behavior to a compact state at a temperature 9 ~ 1.2fo. The collapse temperature 9{x) 
decreases with the asymmetry x = \N^ — A^_|/(A^_|_ + A^_) and vanishes at a; ~ 0.6. For 
fo < 0, a ^-point is present at all x. 
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Polymers in a solvent are subject to both the hard core repulsions between the 
monomers and somewhat longer range (e.g. van der Waals) attractive interactions. At 
high temperatures T, the repulsive interactions are dominant, and the radius of gyra- 
tion (root-mean-squared size) of the polymer scales with the number of monomers A^ as 
Rg ~ N'^ , where u ^ 0.588 in three dimensions. On a discrete lattice, the hard core in- 
teractions are usually modeled by self-avoiding walks (SAWs). The attractive interactions 
are then included by introducing a negative energy for each pair of monomers residing on 
neighboring lattice sites. As temperature is lowered, a point is reached (T = 9) where the 
repulsive and attractive interactions effectively cancel and the polymer behaves essentially 
as an ideal random walk[l][2] with u = 1/2. For T < 9, the polymer collapses into a com- 
pact object with u = 1/3. Numerous Monte Carlo (MC) and exact enumeration studies of 
the ^-point have been performed (see Ref.[3] and references therein). 

The collapse transition for heteropolymers is particularly interesting in view of its 
possible relation to the problem of protein folding[4][5]. While models based on random 
heteropolymers significantly oversimplify the specificity and complexity of real proteins, 
they do bring in fresh perspectives from the statistical mechanics of random systems and 
spin-glasses [6] -[10]. We shall consider polymers formed from two types of monomers, 
labelled by qi = ±1, and subject to a short range interaction 

'^1= 2YI ^^^J^ ^^' ~ ''j) ' (^) 

where r^ are the coordinates of the monomers. Do these interactions modify the collapse 
transition of heteropolymers described above? A perturbative answer is provided by noting 
that [11] that the relevance of such interactions in d space dimensions is controlled by the 
scaling exponent yi = 1 — du. The inhomogeneities are marginal in the compact state 
(z/ = l/(i), and may cause the polymer to choose a particular configuration dependent 
on the specifics of the sequence. The statistical properties of such states are the subject 
of considerable current interest [12] [13]. Here we address the simpler question of how TCj 
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affects tfie non-compact states. For ideal (non-interacting) cfiains witfi v = 1/2, Ti/ is 
relevant in (i < 2: In (i = 1 the polymer swells if like charges repel, and collapses to a point, 
if they attract[ll] (sea also Refs.[14][15]). For all other non-compact states {y > 1/d) weak 
interactions are never relevant. However, we shall demonstrate that strong interactions 
described by Hj do lead to the collapse of a self-avoiding chain. 

We investigate a discretized model in which the only homogeneous interaction is the 
repulsion caused by the constraints of self-avoidance. In the absence of randomness the 
polymer explores all SAWs with equal probabilities. A short range random interaction 
is incorporated by assigning an energy v^qiq^ to every pair of monomers, i and j, on 
neighboring lattice sites. We primarily focussed on vq > 0, for which like charges repel 
and opposite ones attract. This choice corresponds to strongly screened Coulomb inter- 
actions. The asymmetry in the amount of positive and negative charges is measured by 
X = \N^ — N-\/N. We shall show that for moderate values of x the polymer undergoes 
a ^-transition and explore its properties. Negative vq describes a situation in which the 
like monomers attract, and unlike ones repel. A mixture of hydrophobic and hydrophilic 
monomers would exhibit such a tendency. In this case the ^-transition is present for all x. 

We used exact enumeration to study the properties of chains of up to L = 12 steps 
(A^ = 13 monomers). We examined all spatial conformations of SAWs and all possible 
quenched sequences of charges. Taking advantage of rotation and reflection symmetries, 
as well as the degeneracies related to inverting the order of the sequence or the signs of 
all charges, we reduced the number of independent configurations to 4,162,866 for 2,080 
sequences, i.e. a total of roughly 8 x 10^ possibilities. Since the number of cases grows 
by an order of magnitude when L is increased by unity, we could not go beyond L < 12. 
By comparison, the maximal chain-lengths considered in enumeration studies of compact 
walks in d = 3 [13] and d = 2 [16] are about twice longer. In the two dimensional case 
the number of spatial conformations is smaller by 2 to 3 orders of magnitude, while in the 
former case only a random sample of quenches was used. Our results were supplemented by 



MC simulations for L < 95, using the pivoting niethod[17]. (The method provides exceUent 
equihbration at high and intermediate temperatures, but does not permit equihbration of 
almost compact structures[18].) 

Exact enumeration provides a quite detailed picture of the energy landscape. Fig.(l) 
depicts contours of the density of states in the variables i?^ and E for (a) a homogeneous 
chain with Qi = 1 for all i, and (b) a particular quench that is approximately neutral. All 
energy levels are proportional to vq and Fig.(l) corresponds to choosing vq = +1. The 
corresponding densities for i^o < are obtained by simply reflecting the flgures around 
E = 0. From such histograms, the behavior as a function of temperature is deduced as 
follows: At high T the system explores the region of highest density and moves to the 
lowest energy states as T ^ 0. Both chains are SAWs at T ^ oo with R^ ^ 2.4. As 
T is reduced the uniformly charged chain expands; its scaling at T = is still that of a 
SAW, as the only effect of the potential is to exclude conflgurations with chain segments on 
neighboring sites. This effectively increases the range of the excluded volume interaction 
and swells Rg by a flnite factor. The random chain in Fig. (lb) behaves quite differently, 
collapsing to a dense low energy conflguration as T — * 0. 

If the collapse of random chains is through a 9 transition similar to homogeneous 
polymers, we expect that at the transition i?^ ~ L as in ideal random walks. This ex- 
pectation is conflrmed in Fig. (2) which depicts the T-dependence of (quench-averaged) 
R^/L for Ts ranging from 3 to 95. (The quench average is exact for enumeration data, and 
obtained from 20 quenches in the MC simulations.) The curves intersect in the vicinity 
of the same point. Since for T < 6 the polymers are compact {R'^/L ~ L"^'"^), while for 
T > 9 they are expanded SAWs {R'^/L ~ -^'^'^^), the graphs representing larger Ls show 
steeper crossover. Despite a slight L-dependence of the intersection points we can locate 
the ^-temperature with reasonable accuracy at T ^ 1.2^07 but because of small chain 
lengths the crossover exponents cannot be determined reliably. We also examined the heat 
capacity per monomer which exhibits a sharp peak whose height increases with L. The 



position of the peak is close to T = 0.5^0 for short chains and moves towards higher Ts 
for longer chains. 

What happens to the collapse transition if the charges on the chain are not exactly 
balanced (a: > 0)? Due to the excess of repulsive interactions, the chains initially swell on 
lowering temperature. Nevertheless, at sufficiently low temperatures, they may find some 
compact configurations of low energy. The 6 temperature is expected to decrease due to 
charge imbalance. Fig. (3) depicts partial plots of R^/L vs. T at x = 0, 1/2, and 5/6. 
The curves for x = 1/2 are qualitatively similar to x = 0, but with a lower 9 temperature. 
For X = 5/6, however, the curves do not intersect and there is no collapse transition; the 
chain is a SAW at both high and low temperatures. An approximate phase diagram in the 
(XjT) plane is constructed in Fig. (4): the 6'-temperature decreases monotonically with x, 
finally disappearing at a; ~ 0.6. In the T ^ and x -^ 1 limit, the action of minority 
charges is to form loops in the SAW of majority charges. This is somewhat reminiscent of 
a model of SAWs with periodically distributed "stickers" which attach in pairs[19]. Unlike 
our model, the minority stickers are not randomly situated and do not stick to the majority 
monomers. With such stickers, the collapse transition apparently persists all the way to 
X — 1. 

The phase diagram is different for vq < 0, i.e. when like charges attract. The limit 
of X = 1 describes a homogeneous polymer with nearest neighbor attractions. This is 
precisely the model used to investigate the collapse of uniform polymers, and is known 
to undergo a ^-transition at T ^ 3.71^01 [20]. At the other extreme, when x = 0, the 
situation somewhat resembles the case of vq > 0. Although the model is not symmetric 
under changing the of sign of Vq which corresponds to changing the sign of E in Fig. (1), 
the energy landscape of Fig. (lb) appears approximately unchanged by such a refiection. 
Therefore, it is not surprising that there is again a collapse for x = at 9 ~ 1.211)0 1- We 
verified that the model with tJo < does indeed have a collapse temperature for all values 
of X which interpolates between the above two limits. 



Most models of protein folding assume that a homogeneous short range attraction leads 
to formation of a compact phase whose details are then determined by the particulars of 
the sequence. In this work we have shown that even random interactions with no "built- 
in" overall attraction can create a compact state. Our polymers are not long enough to 
determine the critical behavior at the collapse transition. (Even for uniform polymers, 
chains with L > 200 are needed to determine exponents with any degree of accuracy.) 
However, scaling arguments suggest that randomness is irrelevant for ^-chains, and that 
the collapse transition should be in the same universality as that of uniform polymers. 
Some aspects of the model deserve further study. It is plausible that a random polymer 
initially collapses to a dense liquid state, followed by a glass transition upon further cooling. 
(Studies that only enumerate compact states[13] appear to find such a transition.) It would 
also be interesting to investigate the behavior of polymers close to the point at which the 
6'-transition disappears. Detailed investigation of this point, however, requires use of 
much longer chains (L ^ 1000) and methods which permit reliable equilibration of dense 
structures at low temperatures. 

This work was supported by the US-Israel BSF grant No. 92-00026, by the NSF 
through No. DMR-87-19217 (at MIT's CMSE), DMR 91-15491 (at Harvard), and the 
PYI program (MK). YK would like to thank D. Ertas for numerous discussions of the 
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Figure Captions 

Fig. 1: Contour plots of the number of states as a function of R^ (in lattice units) and 
energy (in units of vq) for a 10-step polymer, in (a) uniformly charged case 
(A^_i_ = 11), and (b) for one quench in an almost neutral case {N.^. = 6, A^_ = 5). 
The data bins are of size 1 on the energy scale, and size 1/4 on the R^ scale. The 
smooth contour lines were created by interpolation at levels 0.5, 33, 129, 513, 
2049, 8000, 16000, 32000, 64000. The full circles indicate the location of average 
R^ as the temperature changes between and oo. 

Fig. 2: R^ (in lattice constants) divided by the length of the walk L as a function of tem- 
perature, for A^_|_ = A^_ . Dashed lines represent the results of exact enumeration 
for L = 3, 5, 7, 9, and 11, in order of increasing slope. Full squares and circles, 
open squares and circles, and open triangles represent MC results for L = 9, 13, 
23, 47 and 95, respectively. 

Fig. 3: R^ (in lattice constants) divided by the length of the walk L, for L = 95 (solid 
lines), L = 47 (dot-dashed lines), L = 23 (dotted lines), and L = 11 (dashed 
lines). The groups of lines represent (from right to left) the values of a; = 0, 1/2, 
and 5/6. 

Fig. 4: Phase diagram of a random polymer with t^o > in the plane of temperature (in 
units of vq), and asymmetry x. Vertical bars indicate estimated uncertainties in 
the extrapolated values of the 6'-points. 
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